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ABSTRACT 


A  procedure  is  described  for  determining  dynamic  stress  intensity  factor  histories 
for  a  half  plane  crack  in  an  otherwise  unbounded  elastic  body,  with  the  crack  faces 
subjected  to  tractions  that  result  in  variation  of  the  stress  intensity  factor  along  the 
crack  edge.  The  procedure  is  based  on  integral  transform  methods  and  the  properties 
of  analytic  functions  of  a  complex  variable.  The  procedure  is  illustrated  for  the  case  of  a 
pair  of  opposed  line  loads  suddenly  applied  on  the  crack  faces  along  a  line  perpendicular 
to  the  crack  edge.  An  exact  expression  is  obtained  for  the  resulting  mode  1  stress 
intensity  factor  as  a  function  of  time  for  any  point  along  the  crack  edge.  Some  features 
of  the  solution,  as  well  as  possible  extensions  of  the  procedure,  are  discussed. 
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1.  INTRODUCTION 


In  this  article,  a  procedure  b  described  for  determining  stress  intensity  factor  his¬ 
tories  for  a  class  of  three  dimensional  elastodynamic  crack  problems.  The  geometrical 
configuration  studied  is  a  half  plane  crack  in  an  otherwise  unbounded  solid,  with  the 
crack  faces  subjected  to  tractions  that  result  in  variation  of  the  transient  stress  intensity 
factors  along  the  crack  edge.  Three  dimensional  phenomena  have  received  relatively  lit¬ 
tle  attention  in  dynamic  fracture  mechanics.  The  reflection  of  Rayleigh  waves  travelling 
on  the  face  of  a  crack  and  obliquely  incident  on  the  crack  edge  has  been  studied  (Fre¬ 
und,  1971).  Because  the  fields  are  steady  as  seen  by  an  observer  moving  along  the  crack 
edge  at  an  appropriate  speed,  problems  of  this  type  may  be  reduced  to  two  dimensional 
problems.  This  general  approach  was  also  followed  by  Achenbach  and  Gautesen  (1977) 
in  their  study  of  oblique  reflection  of  plane  harmonic  waves  from  the  edge  of  a  crack. 
Ths  latter  authors  went  on  to  study  harmonic  three  dimensional  loading  of  the  crack 
faces,  obtaining  some  asymptotic  results  for  very  large  or  very  small  frequencies.  In 
some  respects,  the  present  analysis  is  sinoilar  to  the  work  of  Achenbach  and  Gautesen 
(1977).  The  main  differences  are  that  the  present  analysis  allows  for  transient  loading 
and  the  results  are  exact. 

In  section  2,  the  boundary  value  problem  is  stated  for  general  three  dimensional 
crack  face  tractions  that  result  in  the  opening  mode  of  crack  edge  deformation,  or  mode 
I.  The  general  approach  to  be  followed  in  extracting  three  dimensional  stress  intensity 
factors  is  outlined  in  section  3.  In  section  4,  attention  is  restricted  to  the  case  of  a  pair 
of  line  loads  that  are  suddenly  applied  to  the  crack  faces  along  a  line  perpendicv^ar  to 
the  crack  edge,  and  the  results  for  this  case  are  presented  in  section  5.  Some  possible 
extensions  of  the  procedure  are  mentioned  in  section  6. 


2.  FORMULATION 


■:v 

.•i*! 

S’ 


I 


I 

i 


•I 

w' 

■'*V 


I 


Consider  the  elastic  body  containing  a  half  plane  crack  depicted  in  Fig.  1.  The 
body  is  initially  stress  free  and  at  rest.  The  material  is  characterized  by  the  shear 
modulus  ft,  the  Poisson  ratio  u,  and  the  mass  density  p.  The  speeds  of  dilatational, 
shear  and  Rayleigh  waves  are  Cd,  c,  and  Cr,  respectively.  At  a  certain  instant,  tractions 
begin  to  act  on  the  crack  faces,  resulting  in  a  three  dimensional  stress  wave  field  in 
the  material.  The  purpose  is  to  determine  exact  expressions  for  stress  intensity  factors 
as  functions  of  time  and  position  along  the  crack  edge.  To  illustrate  the  approach, 
the  same  transient  normal  pressure  distribution  is  applied  to  each  crack  face,  and  the 
tangential  or  shear  traction  components  are  zero.  For  this  type  of  loading,  the  mode 
of  deformation  is  mode  I  for  each  point  along  the  crack  edge.  Next,  a  mathematical 
formulation  of  a  problem  corresponding  to  the  foregoing  qualitative  description  is  given. 

A  right  handed  rectangular  coordinate  system  is  introduced  in  the  body,  oriented 
so  that  the  z  axis  coincides  with  the  crack  edge,  and  the  half  plane  crack  occupies  y  =  0 
for  z  <  0.  The  normal  traction  on  the  crack  faces  begins  to  act  at  time  t  =  0.  For  the 
time  being,  consider  some  general  traction  variation,  say  Ty(x,z,t)  =  7a-(z,z,t)  on 
y  =  ±0,  where  a-(z,  z,  t)  is  a  given  function  of  position  on  the  faces  and  of  time,  and 
a-  >  0  corresponds  to  tensile  traction.  The  subscript  minus  sign  is  intended  to  imply 
that  the  function  is  defined  for  z  <  0.  The  other  components  of  traction  are  zero,  that 
is,  7’,(z,  z,  /)  =  7’,(z,  z,  t]  =  0. 

In  view  of  the  symmetry  of  the  configuration  and  of  the  applied  traction,  the 
wave  fields  are  expected  to  have  refiective  symmetry  with  respect  to  the  plane  y  =  0, 
and  attention  may  be  limited  to  the  region  y  >  0.  The  displacement  fields  satisfy 
ttx{x,-y,z,t)  =  tt*(z,y,2,f),  uy{x, -y,z,t)  =  -Uy(z,y,2,f),  and  u,(x,-y,  z,t)  = 
tfc(z,  y,  z,  t).  These  properties,  in  turn,  imply  the  symmetry  conditions  (z, 0,  z,  t)  = 


0,  (ryg{x, 0,  Zft)  =  0  and  Uy{x, 0,  z,  t)  =  0  for  a;  >  0  and  for  all  t.  Thus,  the  complete 
set  of  boundary  conditions  to  be  satisfied  by  the  stress  wave  fields  is 


<ryy(x,  0, 2,  t)  =  a-{x,  2,  t)  +  <r+{x,  2,  t) 

<rxy(x,0,2,f)  =0 

(2.1) 

<Ty,(x,0,2,<)  =0 

ttj,(x,0,  2,t)  =  tt-(x,  2,t) 

for  -00  <  X,  2  <  00  and  0<t.  The  definition  of  the  function  <t_,  which  describes  the 
imposed  traction,  is  extended  so  that  (T-  =  0  for  all  x  >  0.  The  function  (t+  represents 
the  unknown  normal  component  of  stress  (Tyy  on  x  >  0  and  =  0  for  x  <  0.  Likewise, 
ti-  is  the  unknown  y  component  of  displacement  of  material  particles  on  the  crack  faces 
for  X  <  0  and  «_  =  0  for  all  x  >  0. 

For  this  system,  the  normal  stress  on  the  plane  y  =  0  ahead  of  the  crack  edge  is 
expected  to  have  the  form 

g+(x,2,<)  =  +0(1)  as  x-^0+  (2.2) 

}/Z7rx 

where  ki(t,  z)  is  the  stress  intensity  factor  history  at  any  point  z  along  the  crack  edge. 
The  main  purpose  here  is  to  determine  this  history. 

The  Helmholtz  representation  of  the  displacement  vector  t?  =  +  V  x  ^  is  intro¬ 

duced,  where  V  is  the  three  dimensional  gradient  operator,  <f>  is  the  scalar  dilatational 
wave  potential,  and  if  is  the  vector  shear  wave  potential.  The  wave  potential  functions 
are  governed  by  the  partial  differential  equations 

c3V2^-^,tt  =  0  c2VV-f.u  =  0  V.^=0  (2.3) 
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in  y  >  0.  The  conditions  that  the  material  is  stress  free  and  at  rest  everywhere  for 
<  <  0  are  expressed  in  terms  of  the  potential  functions  by 

y,  0)  =  ,,  0)  =  =  0  (2.4) 


in  y  >  0.  Finally,  it  is  noted  that  the  boundary  conditions  (2.1)  are  interpreted  as 
though  the  stress  and  displacement  components  were  replaced  by  their  expressions  in 
terms  of  the  potentials  ^  and  if. 


3.  SOLUTION  PROCEDURE 


The  mathematical  problem  described  by  the  differential  equations  (2.3),  the  initial 
conditions  (2.4),  and  the  boundary  conditions  (2.1)  is  linear,  and  integral  transforms, 
along  with  certain  powerful  theorems  from  the  theory  of  analytic  functions  of  a  complex 
variable,  are  used  to  extract  stress  intensity  factor  histories.  In  some  respects,  the 
procedure  is  similar  to  that  introduced  by  deHoop  (1958)  in  his  original  analysis  of  two 
dimensional  systems.  Certain  interpretations  of  the  transformed  fields  at  key  points  in 
the  analysis,  however,  make  it  possible  to  extract  exact  results  for  this  class  of  three 
dimensional  problems.  The  first  step  is  to  apply  the  one  sided  Laplace  transform  on 
time  to  the  differential  equations  (2.3)  and  the  boundary  conditions  (2.1),  taking  into 
account  the  initial  conditions  (2.4).  The  transform  parameter  is  a  and  the  transform 
of  any  function,  say  ^(2,  y,  2,  f),  is  denoted  by  a  superposed  hat,  that  is, 


/oo 


The  parameter  a  is  considered  to  be  a  positive  real  parameter  for  the  time  being. 
Next,  a  two  sided  Laplace  transform  (which  is  equivalent  to  the  more  common  Fourier 
transform  but  which  results  in  slightly  less  cumbersome  algebraic  expressions)  is  applied 
to  suppress  the  dependence  on  z.  The  complex  transform  parameter  is  af,  and  the 
transform  of  any  function,  say  ^{x,  y,  2,  a),  is  denoted  by  the  corresponding  upper  case 
symbol  with  a  superposed  bat, 


/oo 

${x,y,  ^,a)  e-^<^dz 

-00 


Suppose  that  a-  vanishes  for  |2|  >  Zo.  Then  it  is  expected  that 


0(x,  y,  z,s)  =  o  j 


jECt 


as  2  -+  ±00  for  any  c  >  0.  This  condition,  in  turn,  implies  that  the  transform  integral 
in  (3.2)  converges  for  -a  <  Re(f)  <  a  where  a  =  l/cd.  Consequently,  the  integral 
defines  an  analytic  function  in  this  strip  of  the  complex  plane.  From  the  property  of 
uniqueness  of  an  analytic  function  in  its  domain  of  analyticity  (Hille,  1959,  p.  199),  the 
analytic  function  in  the  strip  is  completely  specified  by  the  function  defined  only  on 
the  portion  of  the  real  axis  in  the  strip  -a  <  Re(f)  <  a.  Furthermore,  it  turns  out  to 
be  of  great  advantage  to  restrict  the  range  of  f  to  this  real  interval  in  proceeding  with 
the  analysis.  At  an  appropriate  later  stage,  the  definition  of  the  functions  of  ^  may  be 
continued  into  the  entire  complex  f  plane  by  invoking  the  property  of  uniqueness  of  an 
analytic  function. 

Finally,  the  two  sided  Laplace  transform  suppressing  dependence  on  x  is  ap¬ 
plied.  The  complex  transform  parameter  is  and  the  transform  of  any  function, 
say  4(2,  y,  j,  s),  is  denoted  by  the  same  upper  case  symbol  but  without  the  superposed 
hat,  $(^,  y,  s).  To  be  able  to  continue  with  the  solution  procedure,  some  statement 
on  the  domain  of  convergence  of  the  transform  integrals  over  x  is  required.  If  the  ap¬ 
plied  loading  is  nonzero  for  indefinitely  large  distance  from  the  crack  edge  along  the 
crack  faces,  then  it  may  be  concluded  that  the  integral  over  negative  values  of  x  in  the 
definition  of  ^(^,  y,  f ,  s)  will  converge  for  Re(^)  <  0.  For  x  >  0,  the  wave  fields  are 
zero  beyond  a  cylindrical  wavefront  expanding  from  the  y  axis  at  speed  Cd  with  instan¬ 
taneous  radius  Zo  +  \/x2  +  z^.  To  be  more  specific,  consider  the  elementary  wave  field 
^(x,  y,  2,  t)  =  H{cdt  ~  Zo  -  \/x2  -f  22),  where  H  denotes  the  unit  step  function.  If  the 
Laplace  transform  integrals  over  t,  then  2,  and  then  x  are  formed,  it  is  easily  demon¬ 
strated  that  the  triple  integral  converges  if  (Re(^))*  +  <  a^.  This  condition,  when 

coupled  with  the  convergence  condition  Re(^)  <  0  suggested  above,  provides  the  basis 
for  expecting  that  the  transforms  will  be  analytic  in  the  strip  -^/a^  -  <  Re(^)  <  0 

in  the  complex  (  plane,  with  ^  real  and  in  the  interval  -a  <  f  <  a. 


s 


'I 


It  is  not  possible  to  consider  arbitrary  distributions  of  traction  f)  applied 

to  the  crack  faces.  From  this  point  onward,  attention  is  limited  to  those  distributions 
for  which  the  triple  transform  has  the  form 


Zoo  i-co  *■**’*’  **  (f  >  f ) 


(3.4) 


where  m  is  a  real  number  and  does  not  depend  on  s.  The  condition  (3.4) 

essentially  restricts  a-(x,  x,  t)  to  be  within  a  class  of  certain  separable  or  homogeneous 
functions  of  x,  z,  t. 

The  transforms  are  applied  to  the  partial  differential  equations  (2.3)  to  reduce 
them  to  ordinary  differential  equations  in  y.  The  transformed  boundary  conditions 
(2.1)  are  imposed  on  the  solutions  of  the  ordinary  differential  equations  to  determine 
the  unknown  parameters  of  integration.  As  usual  in  a  problem  of  this  general  type, 
there  are  more  unknown  functions  than  there  are  algebraic  equations  to  determine 
them.  However,  certain  of  the  unknown  functions  are  sectionally  analytic  functions 
(of  ^  in  this  case),  and  the  Wiener- Hopf  factorization  (Noble,  1958,  p.  11)  makes  it 
possible  to  complete  the  solution  by  determining  two  unknown  functions  from  a  single 
equation.  Only  a  few  of  the  intermediate  steps  are  included  here. 

The  solutions  of  the  ordinary  differential  equations  resulting  from  Laplace  trans¬ 
formation  of  the  partial  differential  equations  in  (2.3)  that  are  bounded  as  y  — »  oo 
are 

=  (8.5) 

where  the  parameters  of  integration  P  and  are  unknown  functions  of  their  arguments, 
the  number  m  follows  from  (3.4),  and 

tt  =  <.(?,?)  =  (a«-r-f)'/=,  0  =  =  (3.6) 
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where  a  =  l/cd  and  b  =  1/cg.  The  complex  ^  plane  is  cut  along  f 2  <  |Re(^)|  < 

00,  Im(^)  =  0  so  that  Re(a)  >  0  in  the  entire  cut  ^  plane  for  each  value  of  and 
likewise  for  The  transformation  of  the  condition  V  •  ^  =  0  yields 

^Qx  ~  PQy  d"  fQz  —  0  (3-7) 


where  Qx,Qy,  Qz  are  the  rectangular  components  of  Q.  The  relation  (3.7)  is  obviously 
a  linear  algebraic  equation  for  the  unknown  functions.  The  Laplace  transformation  of 
the  boundary  conditions  provides  four  additional  linear  algebraic  equations  for  the  six 
unknown  functions  F,  Qx,  Qy,  Qz,  U-  and  E+,  where 


j— 00  J—00 

E+(^,f)  =  s"*/  /  a+{x,z,8)e~*^^^'*‘^^Uzdx. 

J—00  J—00 

These  four  equations  are 


(3.8) 


(62  -  -  2f2)F  -  20Qx  +  2 W.  =  (E_  +  E+) 

-2(aP  +  f (Qx  +  ^PQy  +  (62  -  f2  -  2^2)^,  =  0 

2faP  +  (62  -  2f2  -  f eg*  =  0 


(3.9) 


—aP  d-^g®  —  ^Qz  —  U—. 


If  P,  Qx,  Qy,  Qz  are  eliminated  from  the  five  equations  (3.7)  and  (3.9),  the  result  is  a 
single  equation  involving  the  two  remaining  unknown  functions  U-  and  E+,  namely, 


(3.10) 


where  E_(^,  f)  is  given,  and 


f)  =  ‘•(f  +  f ')«(?■  f)y9{?.  f  )  +  (<>'-  2{"  - 


(3.11) 


The  function  defined  in  (3.11)  is  recognized  as  a  modified  form  of  the  Rayleigh  wave 
function.  In  particular,  it  has  zeros  in  the  complex  ^  plane  at  f  =  ±y/c'^  -  where 
c  =  If  Cr  is  defined  by  i?{±c,  0)  =  0,  and  it  has  no  other  zeros  in  the  cut  plane.  Note  that 
(3.11)  is  identical  to  the  standard  Rayleigh  wave  function  when  f  =  0.  Furthermore, 
E+  is  analytic  in  Re(^)  >  -^/a^  -  ^2  and  U-  is  analytic  in  Re(^)  <  0.  The  equation 
(3.10)  holds  in  the  strip  <  Re(^)  <  0.  Thus,  for  any  fixed  value  of  f,  (3.10) 

may  be  solved  by  factorization  in  much  the  same  way  that  an  equation  of  the  standard 
Wiener-Hopf  type  is  solved.  While  the  dependence  on  f  must  be  taken  into  account  in 
the  solution  procedure,  it  requires  no  special  consideration,  and  it  is  simply  understood 
to  be  a  real  parameter  in  the  interval  -a  <  f  <  a  for  the  time  being. 


4.  A  PARTICULAR  TRACTION  DISTRIBUTION 


To  make  further  progress,  it  is  necessary  to  choose  a  particular  applied  traction 
distribution.  A  case  suitable  for  demonstrating  the  general  approach  is  that  of  a  uniform 
line  load  which  suddenly  begins  to  act  along  a  line  perpendicular  to  the  crack  edge. 
Thus,  it  is  assumed  that 

=  (4.1) 

where  H{-)  is  the  unit  step  function  and  £(•)  is  the  Dirac  delta  function.  The  amplitude 
p*  has  physical  dimensions  of  (force) /(length),  and  p*  >  0  corresponds  to  traction  that 
tends  to  separate  the  crack  faces.  Based  on  dimensional  considerations  alone,  it  is 
expected  that  this  loading  will  produce  a  stress  intensity  factor  history  that  depends 
on  z  and  t  according  to 

ki{t,z)  =  p*z-'/^g{tfz)  (4.2) 

Thus,  the  problem  is  to  determine  the  dimensionless  function  5r(')- 

The  ensemble  of  wavefronts  that  results  from  application  of  the  traction  specified 
by  (4.1)  includes  a  variety  of  space  surfaces.  Among  the  wavefronts  are  cylindrical 
dilatational  and  shear  wavefronts  centered  on  the  load  lines  in  x  >  0,  plus  planar 
head  wavefronts  that  intersect  the  dilatational  wavefronts  on  the  surfaces  y  =  ±0  and 
terminate  at  the  line  of  tangency  with  the  cylindrical  shear  wavefronts.  There  are  also 
spherical  dilatational  and  shear  wavefronts  centered  at  the  origin  of  coordinates,  plus 
conical  headwaves,  each  with  its  apex  at  a  point  where  the  spherical  dilatational  wave 
meets  the  z  axis  and  its  terminus  along  a  circle  of  tangency  with  the  spherical  shear 
wavefront.  There  are  also  conical  headwaves  that  intersect  the  spherical  wavefront  on 
the  surfaces  y  =  ±0  for  x  <  0  and  extend  to  circles  of  tangency  with  the  spherical 
shear  wave. 
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If  the  relevant  integral  transforms  indicated  in  (3.4)  are  applied  to  (4.1),  it  is  found 


that  m  =  2  and  that 


The  relationship  between  the  sectionally  analytic  functions  (3.10)  becomes 


where  the  modified  Rayleigh  wave  function  R{^,  g  has  been  replaced  by  the  function 

The  function  5(^,  f)  — ►  1  as  |^1  oo  and  it  has  neither  zeros  nor  poles  in  the  finite  ( 
plane.  The  only  singularities  of  5(^,  f)  are  the  branch  points  at  ^  =  db>/c2  -  ^2  and  at 
(  —  ±\/b^  -  which  are  shared  with  R{(,  f ),  and  it  is  single  valued  in  the  f  plane  cut 
along  <  |Re(^)|  <  ^2-^2^  =  0.  The  property  that  S{1,  f)  =  §(f ,  f) 

for  the  restricted  range  of  ^  may  be  exploited  to  show  that 


where  the  overbar  denotes  complex  conjugate.  The  functions  S+(^,  f )  and  5_($,  f)  are 
analytic  and  nonzero  in  the  half  planes  Re(^)  >  -y/a2  -  fS  and  Re(f)  <  ^a2  -  ^2^ 
respectively. 

With  the  factorization  of  S{i,  j)  complete,  the  relationship  (4.4)  may  be  rewritten 
in  the  form 


P*^+(0,f)  /iKC^_(C,f) 


F±{U)  = 

The  relationship  (4.7)  is  valid  in  the  strip  -y/a^  -  <  Re(^)  <  0.  Furthermore,  each 

side  of  the  equation  is  analytic  in  one  of  the  overlapping  half  planes.  Consequently,  each 
side  of  (4.7)  is  the  unique  analytic  continuation  of  the  other  into  a  complementary  half 
plane,  and  together  the  two  sides  represent  a  single  entire  function  E{0,  parametric 
in  (  (Noble,  1958,  p.  36).  It  remains  to  determine  the  entire  function. 

As  1^1  00  in  the  respective  half  planes,  |F±(^,f)|  =  Furthermore, 

in  view  of  (2.2),  2,  s)  is  expected  to  be  square  root  singular  as  x  0+  for  any 

2,  and  «-(x,  2,s)  is  expected  to  vanish  as  x  — »  0~  for  any  z  to  ensure  continuity  of 
displacement.  Consequently,  from  the  Abel  theorem  concerning  asymptotic  properties 
of  transforms  (Noble,  1958,  p.  36) 

lim  v^S+(^,j)  =  lira  /  d+(x,f,5)e-'f*d2  =  .^/(»,0  (4.8) 

€-»+oo  x-»0+  J—oo 

where 

/®®  »  «  — 

ki  (s,  2)  *  d2 ,  k[Sf  2)  =  lim  \/2ix  a+  (x,  2,  s)  (4.9) 

-00  X-+0+ 

Therefore,  the  left  side  of  (4.7)  vanishes  as  |^|  — ^  00  in  the  left  half  plane,  and  the  same 
conclusion  may  be  drawn  concerning  the  right  side.  According  to  Liouville’s  theorem, 
a  bounded  entire  function  is  a  constant  (Hille,  1959,  p.  204).  In  this  case,  E{()  is 
bounded  in  the  finite  plane  and  E{^)  — »  0  as  |^|  -+  00  so  that  the  constant  must 
have  the  value  zero;  thus,  E{^)  =  0.  This  completes  the  solution  of  the  Wiener-Hopf 
equation  (3.10)  for  the  case  at  hand,  and  the  two  unknown  functions  in  that  equation 


have  been  determined  to  be 


5.  THE  STRESS  INTENSITY  FACTOR  HISTORY  ALONG  THE  EDGE 


It  is  concluded  from  (4.9)  and  (4.11)  that  the  double  transform  of  ki{t,z)y  the 
variation  of  stress  intensity  factor  history  along  the  edge  of  the  crack,  is 


*/(»,?)=?• 


sio;?] 


(5.1) 


where  k  =  2(6^  ~  a®),  as  before. 


Up  to  this  point,  the  variable  f  has  been  understood  to  be  a  real  variable  in 
the  interval  between  -a  and  a.  Because  an  explicit  expression  for  K/(s,()  is  now  in 
hand,  however,  the  definition  can  be  extended  into  the  whole  complex  ^  plane  in  a 
straightforward  way.  The  identity  theorem  for  analytic  functions  (Hille,  1959,  p.  199) 
assures  that  the  extension  or  continuation  is  unique.  When  viewed  as  a  function  of  the 
complex  variable  Ki  is  obviously  analytic  in  the  strip  -a  <  Re(f )  <  a.  Consequently, 
it  may  be  viewed  as  a  convergent  two  sided  Laplace  transform  over  the  spatial  variable 
z.  Furthermore,  extension  of  the  definition  of  f)  into  the  complex  f  plane  reveals 
that  it  is  analytic  everywhere  except  at  branch  points  {•  =  ±a,  ±b,  ±c.  The  function  is 
made  single  valued  by  introducing  the  branch  cuts  extending  from  each  branch  point 
in  a  direction  away  from  the  origin  to  infinity  along  the  real  axis.  The  branches  of 
(o*  -  (ft^  -  and  (c^  -  f2)i/2  ^jtji  positive  real  part  are  selected. 

The  inverse  two  sided  Laplace  transform  of  (5.2)  is 

.  1  rfo+ioo  ^ 

A?i(a,e)  =  5^/  .  (5.2) 

where  (o  is  any  real  number  between  -a  and  a. 

The  central  idea  in  the  next  step  is  to  convert  the  integral  in  (5.1)  to  a  form 
which  will  allow  inversion  of  the  one-sided  Laplace  transform  on  time  by  observation. 


Suppose  that  z  >  0.  Then  the  path  of  integration  is  augmented  by  two  arcs,  quarter 
circles  of  indefinitely  large  radius,  in  the  left  half  plane,  plus  a  path  running  from  the 
remote  point  on  the  negative  real  axis  inward  along  the  upper  side  of  the  branch  cut 
to  =  -a,  around  the  branch  point,  and  then  outward  along  the  lower  side  of  the 
cut  to  the  remote  point;  see  Fig.  2.  These  paths,  together  with  the  inversion  path 
itself,  form  a  closed  contour,  and  the  integrand  of  (5.2)  is  analytic  inside  and  on  this 
contour.  According  to  Cauchy’s  integral  theorem  (Hille,  1959,  p.  163),  the  integral 
of  the  function  along  the  closed  contour  is  zero.  An  immediate  inference  is  that  the 
integral  in  (5.2)  is  equal  to  the  integral  of  the  same  function  taken  along  the  alternate 
path  consisting  of  the  two  arcs  of  large  radius  and  the  branch  line  integral.  According 
to  Jordan’s  lemma,  an  exponential  integral  of  this  type  along  a  remote  circular  arc 
vanishes  as  the  radius  of  the  arc  becomes  indefinitely  large  if  the  function  multiplying 
the  exponential  in  the  integrand  vanishes  uniformly  on  the  arc  as  the  radius  becomes 
indefinitely  large.  In  the  present  case,  |K/(«,f)|  0  uniformly  as  |f|  — >  oo  in  the  left 

half  plane  and,  consequently,  the  value  of  the  integral  along  the  remote  arcs  is  zero. 
Furthermore,  j^/(s,  j)  =  K/(s,  j)  so  that  the  branch  line  integral  may  be  written  as  a 
real  integral, 

^ ^  ^  f  +  jO)}  (5.3) 

A 

where  the  fact  that  Ki{s,  f)  is  to  be  evaluated  on  the  upper  side  of  the  branch  cut  in  the 
left  half  plane  has  been  made  explicit.  In  view  of  the  fact  that  Ki  is  an  even  function 
of  the  integral  takes  the  form 

ki(8,z)  =  -^  Im  I  Ki  (s,  f  +  iO)  |  e"®?  Mf  (5.4) 

and,  indeed,  the  inversion  of  the  one-sided  transform  becomes  obvious.  That  is,  (5.4) 
is  a  product  of  two  transforms,  so  that  Ir/(/,  z)  is  a  convolution  of  the  inverses  of  the 


two  transforms, 


ki{t,z) 


7^  "'•-“I 


(5.5) 


for  ^  >  0,  and  Kj(t,  -z)  =  Kj(t,z).  This  completes  the  process  of  inversion  of  the 
double  transform  for  the  stress  intensity  factor  history  at  each  point  along  the  crack 
edge.  Note  that  the  result  has  the  form  anticipated  in  (4.2).  The  result  is  in  the  form 
of  a  real  integral,  and  some  of  its  features  are  now  examined. 


The  integral  in  (5.5)  cannot  be  evaluated  in  terms  of  elementary  functions.  Its 
general  features  are  nonetheless  evident.  For  example,  the  imaginary  part  of  the  quan¬ 
tity  in  brackets  in  the  integrand  is  negative  for  a  <  (  <  c,  the  it  is  positive  for  f  >  c, 
and  the  integral  is  singular  at  t/2  =  e.  The  integral  has  been  evaluated  numerically, 
and  the  result  of  numerical  evaluation  for  a  Poisson  ratio  of  i/  =  0.3  is  shown  in  Fig.  3. 
Knowing  the  asymptotic  behavior  of  the  generalized  Rayleigh  wave  function,  namely, 
that  iJ(0,  f )  -*  as  j  -*■  -f  00,  the  integral  (5.5)  implies  that 


Urn  k/(f,  z)  =  k/(oo,  z)  = 

i—*00 

I 


JL 

y/Wz 


(5.6) 


for  any  fixed  2:  >  0.  The  result  of  the  limiting  process  is  consistent  with  the  indepen¬ 
dently  calculated  equilibrium  stress  intensity  factor  distribution  along  the  crack  edge 
for  this  configuration  (Tada  et  al,  1973,  p.  23.1).  The  stress  intensity  factor  history  in 
Fig.  3  for  any  fixed  value  of  z  is  normalized  by  its  long  time  limiting  value  k/(oo,  z). 

The  general  features  of  ki{t,z)  derive  from  the  nature  of  the  wave  fields.  For 
any  value  of  2  >  0,  the  stress  intensity  factor  is  zero  up  until  the  arrival  of  the  first 
dilatational  wave  at  time  t  =  zjcd.  Upon  sudden  application  of  the  compressive  line 
loads  on  the  faces  of  the  crack,  the  initial  response  is  dilatational  and  the  surfaces  at 
points  adjacent  to  the  line  loads  tend  to  bulge  outward.  The  first  wave  arriving  at  the 


observation  point  along  the  crack  edge  carries  this  outward  bulge.  The  crack  faces  tend 
to  move  toward  each  other,  and  this  tendency  is  reflected  in  the  tendency  for  the  stress 
intensity  factor  to  become  negative  following  arrival  of  the  leading  dilatational  wave. 
This  effect,  which  is  similar  to  that  discussed  by  Freund  (1974)  in  a  two  dimensional 
situation,  persists  until  the  arrival  of  the  Rayleigh  wave  at  time  t  =  zfcr  (cj/cr  =  2.02 
for  V  =  0.3).  The  stress  intensity  factor  history  is  logarithmically  singular  at  this 
instant,  and  it  tends  to  increase  thereafter.  The  transient  stress  intensity  factor  history 
kj{ty  z)  approaches  its  long  time  limit  kj  (oo,  z)  quite  slowly,  with  the  ratio  of  the  former 
to  the  latter  being  about  0.81  at  Cdifz  =  20  and  about  0.92  at  Cdtfz  =  100. 


6.  CONCLUDING  REMARKS 


This  completes  the  analysis  of  the  three  dimensional  stress  intensity  factor  history 
for  the  particular  applied  traction  distribution  (4.1).  Results  could  be  derived  for  a 
number  of  other  traction  distributions  in  a  similar  way.  For  example,  suppose  that 
crack  face  traction  in  the  form  of  a  pair  of  opposed  concentrated  forces  of  intensity 
p*,  one  acting  on  each  face  of  the  crack,  begins  at  the  crack  edge  and  then  moves  at 
uniform  speed  v  in  the  negative  x  direction.  In  this  case,  (4.1)  is  replaced  by 

=  -p*6{x  +  vt)6{z)H{t)  (6.1) 

Similarly,  if  a  pair  of  line  loads  of  intensity  p*  acts  over  x  =  0,  -vt  <  x  <  0  so  that  the 
step  moves  at  uniform  speed  v  in  the  negative  x  direction, 

ff_(x,x,<)  =  -p*H{x  +  vt)6{z)H{t)  (6.2) 

Yet  another,  slightly  more  complex,  illustration  is  obtained  if  a  pair  of  line  loads  of 
intensity  p*  begins  to  act  along  lines  inclined  at  an  angle  <jj  to  the  x  direction.  In  this 
case, 

a-(x,z,t)  =  -p*5(xsina;  + ^cosa;)£f(-x)/f(t)  (6.3) 

I  and  so  on.  The  present  analysis  opens  the  way  for  examination  of  these  and  other 

boundary  value  problems  connected  with  three  dimensional  phenomena. 

I  A  far  more  challenging  situation  is  presented  if  an  opposed  pair  of  concentrated 

forces  of  intensity  p*  suddenly  begins  to  »:t  on  the  crack  faces  at  a  point  at  some  fixed 
distance  (>0  from  the  crack  edge,  that  is, 

I 

1 


a-(x,z,t)  =  -p*6{x  +  C)6{z)H{t) 


(6.4) 


The  means  of  obtaining  an  exact  stress  intensity  factor  solution  in  this  case  has  been 
elusive.  The  equivalent  two  dimensional  problem  was  solved  by  means  of  the  dynamic 
weight  function  method  (Freund  and  Rice,  1974).  In  principle,  the  weight  function 
method  is  also  applicable  for  this  three  dimensional  situation.  However,  the  procedure 
requires  that  the  dependence  of  tt(x,  z,  s)  on  x  and  z  be  known  quite  explicitly,  and  the 
way  to  extract  this  dependence  from  the  double  inversion  integral  (4.11)2  has  not  yet 
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FIGURE  CAPTIONS 


Fig.  1  Geometrical  configuration  of  the  elastic  solid. 

Fig.  2  The  complex  f  plane  showing  the  singularities  of  ki{s,  ()  and  the  integration  path 
for  evaluation  of  (5.2). 

Fig.  3  Normalized  stress  intensity  factor  ki{t,z),y¥zfp*  versus  Cdtfz. 


Figure  2 


